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\ Abstract 



We derive the complete set of velocity-enhanced terms in the expansion of the total cross sec- 
tion for heavy-quark pair production in hadronic collisions at next-to-next-to-leading order. Our 
expression takes into account the effects of soft-gluon emission as well as that of potential-gluon 
exchanges. We prove that there are no enhancements due to subleading soft-gluon couplings 
multiplying the leading Coulomb singularity. 



Oh 



1. Introduction 



Hadronic heavy-quark pair production is of interest not only because of its phenomenological 
\Q \ relevance in the particular case of top quarks, but also because of the theoretical insights that 

may be gained into the singularity structure of QCD in the presence of massive partons. 

A complete result for the next-to-next-to-leading order (NNLO) 0(af) corrections is still elu- 
sive despite substantial recent progress. First, the virtual corrections at NNLO were obtained 
in the limit of large invariants s, —i Am 2 0, Q thanks to the understanding of this limit as 
an alternative regularization scheme for collinear divergences 0, 0]. Subsequently, the complete 
amplitude in the case of quark annihilation was derived with numerical methods [5|, followed by 
partial analytic results @ 0]. The latest studies of the structure of massive gauge amplitudes 
[1, E3, El, E EH led to the derivation of the soft and collinear divergences of the gg(qq) — > QQ 
two- loop amplitudes in analytical form [Hj |. Despite all this, the NNLO program for the calcula- 
tion of the heavy-particle pair production cross section will only be com plet ed once these two-loop 



corrections have been combined with squared one-loop corrections [14 , Il5l . Il6j , the one- loop cor 



rections to tt + j [17], and the real corrections. 



In this note we build on our recent work on soft-gluon radiation [1 21 ] and soft-gluon radiation 



in the presence of Coulomb enhancements [ll| to derive the complete set of velocity-enhanced 
terms in the expansion of the total heavy-quark pair production cross section at NNLO. The 
result contains terms of the form f} 1 log 5 /3, with j3 = y/X — 4m 2 /s the velocity of the heavy quark, 
and —2 < i < 0, < j < 2i + 4, which we provide, apart from the constant term i = j = 0, 
which remains inaccessible with the methods used. As discussed in [ll| , the In (3 coefficient at 
NNLO receives contributions from heavy-quark potentials other than the Coulomb potential, and 
from soft-gluon effects not contained in the standard resummation formula for the total cross 
section due to subleading soft-gluon couplings multiplying the leading Coulomb singularity. We 
compute these terms below, and prove that the subleading soft-gluon effects vanish for the total 
cross section. Our results may be used for improving approximate NNLO top-quark production 
cross section calculations [HI . We also provide a general formula for the velocity-enhanced terms 
at NNLO for arbitrary colour representations of the particles involved, which requires as only 
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process-specific input the constant term of the NLO virtual amplitude at threshold in every colour 
and spin channel. 



2. Sources of enhancement 

There are two general sources of enhancement of the partonic cross section near threshold. One 
is connected with the emission of soft gluons, resulting in up to two powers of ln/3 per emission, 
whereas the other is due to potential exchanges of gluons (the Coulomb potential being the most 
prominent example), yielding up to one factor of 1/(3 and ln/3 per loop. 

The NNLL soft-gluon effects can be readily derived either by inverting the Mellin-space trans- 
form of Ref. [12j or, directly in x-space, from the resummation formula of Ref. [ll|. The latter 



approach has been studied before [lfj, l20| in the context of the Drell-Yan and Higgs-boson produc- 
tion processes at hadron colliders. The soft-gluon enhancement of the Coulomb-potential effects 
has been considered at NLO/NLL in Ref. [21( assuming factorization of the two effects. A priori, 
there could be a highly non-trivial intcrtwinement between the two. T he g eneral factorization of 
leading soft-gluon and Coulomb effects has been studied recently in Ref. [llj to all orders, resulting 
in the formula 
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for the partonic cross sections, which justifies the multiplicative factorization of the two effects in 
Mellin space to NNLO. In the case at hand, where we are only interested in the NNLO expansion 
of the resummed result, we can describe the problem as a sum of a pure NNLO soft-gluon exchange 
(which requires the knowledge of the two- loop anomalous dimension 11, E2]) contained in the two- 
loop contribution to the soft function W^, a ; a soft-gluon enhancement of the Coulomb-potential 
exchange due to the convolution of the one-loop terms in Jr o and W^, a [22| : and, finally, two-loop 
non-Coulomb potential and kinetic-energy corrections as described by non-relativistic effective 
theory (NRQCD), which for the present purpose may be thought of as two-loop contributions to 
the non-relativistic function J Ra . By factorizing the "hard" cross sections Hai (m, /x) from the soft 
and Coulomb effects, one generates not only logarithms of the type 1//3 x log 2 (3 and 1/(3 x log/3, 
but also a process-dependent non- logarithmic term proportional to 1/(3, due to the product of the 
matchin g co efficients and the Coulomb potential. Our result below differs from the one given in 
Ref. 18, 2311 due to this effect, the value of the two-loop soft anomalous dimension 11, 12j], and 



231 ] due to this effect, the value of the two-loop soft anomalous dimension [11. 
the ln/3 terms from the non-Coulomb effects, which we now derive. 

We present two ways to obtain the desired result. The first is based on an explicit calculation 
of the potential contribution in NRQCD. We first generalize the expression for the colour-singlet 
heavy-quark potential in momentum space given in Ref. 24] to arbitrary colour representations 
R a . To obtain the NNLO ln/3 terms it is sufficient to use the four-dimensional potentials. The 
required terms read 
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where D Ra is the strength of the Coulomb potential in representation R a , such that D Ra = —Cp 
for the singlet and D Ra = —(Cp — Ca/2) for the octet representation (note the sign convention!), 
and w S pin = and —2/3 for a it pair in a spin-singlet and spin-triplet state, respectively. The 
non-Coulomb potentials, including a new result for the one-loop l/r 2 -potential in an arbitrary 
representation, are those in the second line. The first line refers to the Coulomb potential and 
its one-loop correction that is already dealt with as described above. To NNLO we need the 
resummed insertions of the non-Coulomb potentials used in the e + e~ — 1 



tt calculation of Ref. [24| , 
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Figure 1: Graphs relevant to the 1/r 2 potential contributions in the singlet and octet channels discussed in the 
text. Crosses correspond to the singlet or octet colour-projection operators. 



given explicitly in Ref. [25| . expanded to NNLO, which results in very simple expressions. Including 
the relativistic kinetic-energy correction we find the non-Coulomb contribution 

ax\nc = «'(M 2 ) ln/3 [-2Dl a (1 + %>in) + D Ra C A ] (3) 

to the total cross section, with the Born cross section in the spin and colour channel X. For 
top quarks the Born cross section in the qq initiated channel is a pure colour-octet spin-triplet, 
whereas in gluon-gluon fusion the tt state is spin-singlet but colour-octet or -singlet. 

The second derivation of the non-Coulomb logarithms uses known results on the threshold 



expansion of the e + e — > tt [26| and 77 — > tt 2jj processes at NNLO. The potential contributions 



are implicit in these results, and the two processes cover the singlet-triplet spin dependence of the 
results for the enhanced terms in exact correspondence to the hadronic case. The only non-trivial 
issue is the colour dependence since we also need the colour-octet case. It is well known that for 
an interaction with colour structure T a ®T a , the transition between singlet and octet is obtained 
by a simple change Cf — > Cp — Ca/2 corresponding to the different value of Dn a . But the 1/r 2 
potential comes also from exchanges of two gluons as depicted in Fig. [1] An explicit check proves 
that each of the diagrams gives the correct contribution (as far as colour is concerned), with the 
same replacement as before. Thus, we obtain the correct results for hadronic tt production by 
keeping only the velocity-enhanced terms from the respective formulae of [26l . I271 ]. making the 
replacement Cf —* Cf — Ca/2 for colour-octet contributions, and removing the contribution 
from the hard matching coefficient at one-loop multiplying the one-loop Coulomb potential. The 
latter step is crucial in obtaining the correct result, since the appropriate matching coefficients 
corresponding to the processes considered have already been taken into account in the soft gluon 
enhancement of the Coulomb contribution as described above. 

There could be another enhanced single or double logarithm of velocity at NNLO from the 
product of a a s / (3 Coulomb term multiplying an a s /31n 2 /3 or a s (3 In (3 term from a beta-suppressed 
subleading soft-gluon coupling [ll[ . Such suppressed couplings exist for the emission of soft gluons 
from the initial state as well as from the final state. We now show that such terms do not appear 
in the total pair production cross sectionQ To this end we imagine obtaining the cross section 
by evaluating the imaginary part of forward-scattering graphs such as those of Fig. O We first 
consider the subleading coupling to the heavy-quark loop, so the gluon coupling to the external line 
in Fig. [5] is the standard eikonal coupling. In the framework of non-relativistic effective theory the 



subleading gluon coupling corresponds to the x ■ E interaction [29|, |30( . An expansion of the heavy- 
quark loop in the velocity can be extracted directly by the strategy of regions [3l| . As described 
in the latter work, it is sufficient to consider the following regions of integration momenta in the 
partonic cms frame where the sum of the heavy-particle momenta is (2m, 0): hard [k ~ to, with to 
the heavy-quark mass and k a loop momentum), potential (k ~ mf3 2 , k ~ to/3), soft (fc ~ to/3, 
k ~ m/3) and ultrasoft (k° ~ to/3 2 , k ~ to/3 2 ). At NNLO the diagrams corresponding to the soft 
region vanish, as they generate only scaleless integrals. The source of singular terms is in the 
potential and ultrasoft regions. By the velocity scaling only terms corresponding to the potential 
three-momentum contribute odd powers of /3. The same scaling arguments also show that, in any 



1 This can be anticipated from the known expansion of the NLO cross section 28]. From this result one can 
readily verify that the logarithms of velocity appear in terms suppressed by even powers of j3 relative to the leading 
terms, i.e. no a s /31og/3 terms of the mentioned origin are generated. 



3 



Figure 2: Example graphs with contributions from both the ultrasoft (gluons depicted with wavy lines) and potential 
(gluons depicted with dashed lines) regions. The crosses denote effective interactions, the structure of which is 
irrelevant to the argument of the text. 

denominator containing a combination of a potential and an ultrasoft momentum, the ultrasoft 
momentum will be (multipole) expanded. Therefore, the denominators containing potential three- 
momenta will not depend on the direction of any external three-momentum (unlike denominators 
containing an ultrasoft three-momentum). In consequence, rotational invariance implies that all 
integrals with an odd number of potential three-momenta in the numerator vanish. Thus, given 
a term with a specified power of /?, the next higher-order contribution will be suppressed by a 
relative factor of (3 2 , smaller than the terms we seek. 

Next, regarding the subleading soft-gluon couplings to the initial state, the relevant expan- 
sion is one in transverse momentum. The effective Lagrangian for the corrections to the eikonal 
app roximation is given in soft-collinear effective theory by £ (x^_n v _ W c gF^W^) £ for quarks 



32 . |33j , and similar terms involving transverse derivatives or factors of x± for the couplings to 
collinear gluons, and of soft quarks. None of these terms can contribute a beta-suppressed term, 
since the initial-state momenta in Fig. [2] can always be chosen to have zero transverse momen- 
tum, implying that loop integrals with transverse-momentum factors in the numerator vanish by 
arguments similar to those applied to the heavy-quark couplings. This completes the proof, that 
we have correctly taken into account all possible sources of singular terms in the expansion of the 
cross sections for heavy-quark pair production at NNLO by including the extra terms from the 
non-Coulomb potentials. 

Note that some of the cuts of Fig. [2] correspond to three-particle colour correlations at the 
amplitude level, for which the infrared divergence structure has recently been given in Ref. (l3| . 
The latter work shows that the infrared-singular three-particle correlations may not vanish in 
the limit (3 — > in the amplitude, but that they do in the virtual contributions to the total 
cross section at NNLO in the particular case of top quarks because of colour projections 12, 13| . 



Our arguments above prove that there are no contributions to the In terms from three-particle 
correlations in both, the virtual and real corrections. This holds independent of particular colour 
representations for purely kinematic reasons. 



3. Results 



Next we present the main result of this paper, namely the expansion of the two-loop partonic cross 
section close to the partonic threshold (3 = 0. As we emphasized above, our result is complete up to 
the so-called constant termf0 C^q. 



that goes beyond the scope of the present work. Setting hr 
cross-section close to threshold reads: 
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2 This standard terminology is somewhat misleading in this process. Due to the non-trivial (3 dependence of the 
Born cross section, the contribution of the "constant" term to the cross section is, in fact, proportional to /3. 
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where I = 1,8 is a colour index and ij = {qq,gg), whereas a s (/u 2 ) is defined in the MS scheme 
with ni (number of massless quarks) flavours. The derivation of the coefficients of ln™(/Lt 2 /m 2 ) 
with n — 1,2 from one- loop results and splitting functions is given in | Appendix B| The non-trivial 
scale-independent two-loop contributions <J^'^ read: 
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In order to construct the two-loop correction to the colour-averaged cross section from the 



colour-state specific components given above, one has to first multiply the two-loop contributions 
Eqs. (|5l7p and ^ by, respectively, <j^\ and cr-^ (see Eq. Q), and then add them together. The 
singlet/octet Born terms can be found in | Appendix B[ Finally, by setting = m, all colour factors 
to their numerical values, and rii = 5 as applicable to top-quark production, we get the following 
result for the colour-averaged total inclusive cross-section close to partonic threshold: 
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which differs in the coefficients of the l/f3 and ln/3 terms from the expressions given in [23| for 
the reasons mentioned in section [51 

In conclusion, the above formulae contain all velocity-enhanced terms in the total hadronic pro- 
duction of heavy quarks at NNLO near the partonic threshold. A compact general result for the 
velocity-enhanced terms in the production of equal-mass heavy-particle pairs in the collisions of 
massless particles for arbitrary colour representations is provided in appendix A. 
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Appendix A. A formula for arbitrary representations 

Here we provide the velocity-enhanced terms at NNLO for the production of a pair of heavy 
particles with equal mass m in the scattering of massless partons in colour representations r 
and r', respectively, under the assumption that the Born cross section (which is factored out 
as in Eq. 0) admits an S'-wave term proportional to (3. The heavy-particle pair is in colour 
representation R a and a definite spin state. The threshold expansion reads 
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We obtained this result by expanding the resummation formula |T]) to NNLO and adding the non- 
Coulomb terms according to Eq. ([3]). C' r , C T i and C Ra denote the quadratic Casimir operators 
of the colour representations, /?o = ^-Ca — ^niTf is the one-loop beta-function coefficient, and 
Ls = ln(8m//x). The quantities D Ra , a\ = ^-Ca ~^§~ n iTf and v sp - ln are connected with the heavy- 
quark potentials as discussed in the main text. is the one-loop hard matchin g co efficient in 
the resummation formula for hadronic heavy-particle pair production at threshold , when the 
heavy-particle pair is in colour and spin state X. Alternatively, it can be deduced from the 
constant term in the threshold limit of the NLO production cross section in colour and spin 
channel X by comparing the expansion of to the formula 

2n 2 D Ra 



7 (1) 







+ 4(C r + CV') 



[m 2 ( 8w/32 ) 


llTT 2 ! 

+ 8 ~ 24 J 







■i(C Ra +4(C r + C r ,)) In 



8mf3 2 



12C Ra +2Re[C^ 1) ]+0(/?) 



(A.2) 



The results for tt production in the main text can be generated from the general formula by 
inserting the relevant colour and spin factors. The required matching coefficients ReC^ can be 
determined from [34, 35 1 (for convenience of the reader, the NLO cross sections are reproduced in 
Eq. (|B.6[) below) and (|A.2[) . Note that the cubic ln/x dependence in the L| term in (|A.1[) cancels 
with a corresponding term in the product Re [C^p] L s as required since can depend on ln/i 
at most quadratically. 



Appendix B. Derivation of the scale dependence 



The scale-dependent terms in Eq. (J4J) can be obtained from Eq. (|A.ip . Here we give an independent 
derivation from the known LO and NLO results; the procedure has been detailed, for example, 
in Ref. [36| and we adopt it in the following. To the best of our knowledge, these terms are not 
available in the literature in completely analytical form. 

We consider the case /i R = fip = fi. To simplify the following discussion, we introduce the 
functions s 
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Note that as follows from Eq. (|B.5p below, the overall factor £t].j/j9 is a beta independent constant. 
It is factored out for later convenience. We will not consider the ij ~ qg subprocess since cr qg j = 
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0((3 3 ). The scaling functions appearing in Eq. (|B.1[) read: 
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In the above equation we need to consider the cross sections as functions of the dimensionless 
variable p = 4m 2 /s, and with = y/1 — p. We have also introduced the "hat" notation such 
that for any function a{p) we have &{p) = a(p)/p; for its origin see the discussion in Ref. [36}. 
Also, summation over repeated indexes k = q, q, g is understood and 0o and 0\ are the QCD beta 
function coefficients 

00 = y0 A - jp, , 01 = y C\ - yC #i - 2C Fni , (B.3) 

with C A = 3 and C F — 4/3. Throughout this letter we take the number of active quarks to be 
equal to the number of light quarks n,. The functions are the one- and two- loop space-like 
DGLAP evolution kernels defined as an expansion in a s /(2n). 

For the reactions of interest, and neglecting terms that contribute beyond O(0), Eqs. (|B.2j) 
reduce to: 



= Po&( P )-2pW®sf£(p), 
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The results for the gg initiated subprocess can be obtained from Eqs. (|B.4[) by replacing everywhere 



the pairs qq and qq by gg, see also Ref. [37 1 



The above convolutions can be performed in a straightforward manner, by expanding the 
splitting functions around x — 1 and keeping only the 1/[1 — terms and 8- functions. To make 
our presentation self contained, in the following we present all terms from Eq. Q that are not 
already given in Eqs. (|5I6I7[) . We do not present the terms 0(1), i.e. the terms that are not 
enhanced by inverse power of or by a power of In 0, since they are at the same level as the 

(2) (21 (2) _______ 

presently unknown coefficients Cqq ,Cg g x and C gg \a appearing in Eqs. (|5l6l7p . Nevertheless, we 
have calculated these 0(1) terms and they can be found in electronic form available with the 
preprint of this paper. 

The leading-order (Born) terms appearing in Eq. ((4]) read: 
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with TV = 3 the number of colours. The NLO terms appearing in Eq. ((J) read: 

<r£',8 = f C A -|n z +C F (10- 161n2)-16C F ln/3, 

2 

"X.b = (20 F -C j4 )^ + 32C F ln 2 /3+[32C F (-2 + 31n2)-8C A ]ln/3 
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T (i.o) _ 



r (i,i) _ 

T (l,0) _ 



-^C F (-24 + 7T 2 + 63 In 2 - 48 In 2 2) + ^C A (308 - 9vr 2 - 180 In 2) 

4 32 
+-n,(-5 + 61n2)-— , 

-16C A (-l + ln2) - 16C A ]n/3, 

7T 2 

(2C F - C^) — + 32C A In 2 (3 + 2AC A (-3 + 4 In 2) In /? 

P 

+C F {-2Q + n 2 ) + ^(7^(504- 17tt 2 - 624 In 2 + 384 In 2 2) , 

(1,1) 



a 



2C F ^- + 32C A In 2 + 32C A (-2 + 3 In 2) ln/3 

+C* f (-20 + tt 2 ) + ^6^(204- 7tt 2 - 288 In 2 + 192 In 2 2) . 
3 

Finally, the NNLO terms in Eq. (|4]) proportional to ln™(^ 2 /m 2 ) read: 



(B.6) 



T (2,2) 
T (2,l) 



T (2,2) 
T (2,l) 



T (2,l) 



ln/3 + 0(l) : 



128C F 2 ln 2 /3 + 32C F 2 (-5 + 81n2)- ^C F C A + ^-C^ 

7T 2 

— [(WC A C F - 32C F 2 )ln/3 - 11C/ - 12C F 2 - 4C F n, + C A (28C F + 2n,)] 

-512C F 2 In 3 (3 + [480C F C A - 64C F n z - 64C F 2 (-21 + 32 In 2)] In 2 [3 

64 

-88C A 2 + 16C A n, + — C F 2 (-102 + 7tt 2 + 156 In 2 - 120 In 2 2) 



/ 512 32 fS'Jj \ 

+c F ( ^- -^(- 17 + 22 ln2 ) + -tCaC- 136 + 27r2 + 141 ln2 ) ) 

9 3 3 / 



32 



ln/3+0(l), 



128C/ In 2 /3 + 

.2 



8 16 

-C A 2 (-107 + 961n2) + y C A n i 



ln/3 + C(l). 



7T 



(16C A 2 - 32C A C F )ln/3 - ^-C A 2 + C F (^-C A - ~n ; J + ^C A n t 



-512C/ \n 3 (3 



32 



64 



+ 



- y C A 2 (-167+ 1921n2)- y C A n ; 
16 



In 2 /? 



-16C F C A (-20 + 7r 2 ) - —C A n; (-37 + 36 In 2) 



-C/(-4391 + 201tt 2 + 5292 ln 2 - 2880 In 2 2) 
9 



ln/3 + 0(l) 



r (2,2) _ (2,2) 



7T 

1 I 



22 



-32C F C A kip+—C F Cj 



:C 



F'H 



512C A 2 In 3 (3 



64 



32 „ 2 



C/(155- 192 ln 2) 



ln 2 /3 



16C F C A (20-7r 2 



32 



C A n;(17- 181n2) 



16 



C A 2 (-1963 + 96tt 2 + 2502 ln 2 - 1440 In 2 2) 



ln/3 + 0(l). 



(B.7) 



Eqs. (|B.7[) are in agreement with the partially numerical results provided in Ref. [2^ 
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